A geometric proof of the equality between entanglement and edge spectra 
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The bulk-edge correspondence for topological quantum liquids states that the spectrum of the reduced density 
matrix of a large subregion reproduces the thermal spectrum of a physical edge. This correspondence suggests 
an intricate connection between ground state entanglement and physical edge dynamics. We give a simple 
geometric proof of the bulk-edge correspondence for a wide variety of physical systems. Our unified proof 
relies on geometric techniques available in Lorentz invariant and conformally invariant quantum field theories. 
These methods were originally developed in part to understand the physics of black holes, and we now apply 
them to determine the local structure of entanglement in quantum many-body systems. 



INTRODUCTION 

Since the discovery of the fractional quantum Hall fluids 1 
and the subsequent elucidation of their topological structure 
2, 3, it has become clear that entanglement plays a crucial role 
in wide variety of zero temperature quantum phases of matter. 
To describe such phases of matter, it is necessary to under- 
stand their pattern of long range entanglement since they fail 
to be distinguished by any symmetry breaking pattern. The 
basic measure of entanglement in the ground state is the entan- 
glement entropy of a subsystem defined as the von Neumann 
entropy of the reduced density matrix of the subsystem. The 
boundary law for entanglement entropy states that the entropy 
of a subsystem of linear size L in the ground state typically 
scales like i.e. like the boundary of the region (d is the 
space dimension). Ref. 4 provides a comprehensive review of 
this fundamental result. Entanglement entropy has since been 
used in Refs. 5, 6 to give meaning to the notion of long range 
entanglement in topological fluids. Entanglement considera- 
tions have also led to a revolution in our understanding of Id 
physics 7-9, promising variational states in higher dimensions 
10-12, and a classification of Id phases 13, 14. At present, 
we are still continuing to explore the role of entanglement in 
quantum many-body physics. 

More recently, Li and Haldane in Ref. 15 drew attention to 
interesting physical information encoded in the full spectrum 
of the reduced density matrix of which the entanglement en- 
tropy is a single measure. Of course, the full spectrum is com- 
plicated and for a large enough subsystem, impossible even 
to list effectively, but what Li and Haldane argued was that 
the entanglement spectrum contained a universal part, even at 
moderate size, that was characteristic of the phase of inter- 
est. To understand their claim, we write the reduced density 
matrix of region R as pR = e R , where we have defined 
the "entanglement hamiltonian" Hr. This statement by itself 
contains very little information. The boundary law implies 
that Hr behaves, at a very crude level, like the Hamiltonian 
of a lower dimensional system. What Li and Haldane argued 
was that for quantum hall systems the universal part of Hr 
was actually given by the true dynamical Hamiltonian of a 
physical edge. Thus the entanglement cut becomes a physical 
cut, and the ground state informs us about quantum dynam- 



ics at the edge. The entanglement spectrum has since been 
studied intensively, see Refs. 16-18. We should also mention 
that information equivalent to the Li-Haldane proposal can be 
extracted from generalizations of the von Neumann entropy 
called Renyi entropies which have been studied previously. 
These entropies are defined as S a (R) = j^- In (tr(p Q )) and 
can be thought of as computing the entanglement entropy as a 
function of "entanglement temperature" 1/a. 

Since the Li-Haldane proposal, known as a bulk-edge cor- 
respondence, Refs. 19-23 have offered a variety of proofs of 
the conjecture and variations of it in various situations ranging 
from non-interacting topological insulators to fractional quan- 
tum Hall fluids in 2 + 1 dimensions. There were even hints 
of a bulk-edge correspondence in the early work of Ref. 5. 
In this paper, we present a proof of the bulk-edge correspon- 
dence for a wide variety of physical systems. Our proof has 
the advantage of simplicity, physical transparency, and gener- 
ality. It is also appealingly geometrical in nature. Our main 
technical tools are a set of powerful geometrical constructions 
in Lorentz invariant and conformally invariant field theories 
that relate the entanglement spectrum of special subsystems 
to thermal spectra in appropriate spacetimes. 

We make use of existing techniques in conformal field the- 
ory described in Refs. 24-26 to rigorously establish the in- 
tuition that the boundary used to define entanglement entropy 
functions much like a real boundary in topological liquids. In 
particular, if at a real boundary the system has gapless edge 
modes, then we argue that the entanglement spectrum also 
contains signatures of these gapless modes. Thus we estab- 
lish very generally a powerful link between entanglement and 
quantum dynamics, connecting ground state properties with 
edge dynamics. Our results treat systems in a variety of di- 
mensions with and without interactions in a completely uni- 
fied framework. We provide a new derivation of topological 
entanglement entropy in 2 + 1 dimensions, a proof of the bulk- 
edge correspondence for entanglement spectra in fractional 
quantum Hall states, a proof of the bulk-edge correspondence 
for topological insulators in 3 + 1 dimensions, and a proof of 
the bulk-edge correspondence for fractionalized topological 
insulators in 3 + 1 dimensions. 

This paper is organized as follows. The first section con- 
tains an introduction and motivation. Section two describes 
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our geometric tools and the intuition behind them. Section 
three treats the original Li-Haldane situation of fractional 
quantum Hall fluids followed by analogous calculations for 
topological band insulators and a class of recently proposed 
fractional topological insulators. Finally, we conclude with 
some comments on future directions. 



ENTANGLEMENT HAMILTONIANS FROM GEOMETRIC 
FLOWS 

Our goal is to understand the spectrum of the reduced den- 
sity matrix of a disk (or more generally, a ball) in various kinds 
of topological quantum liquids. Our main tool is a result from 
relativistic conformal field theory that identifies the spectrum 
of the "entanglement hamiltonian" or "modular hamiltonian" 
(as it is known in the field theory literature) of a ball in d 
spatial dimensions with the spectrum of the usual dynami- 
cal hamiltonian of the conformal field theory on a hyperbolic 
spacetime of the form H d x R. As a geometric aside, for 
d = 2, H 2 is the familiar Poincare disk which can be covered 
by coordinates x, y with x 2 + y 2 < 1 and has the metric 



ds 2 = 



dx 2 + dy 2 



(1 - (x 2 +y 2 )) 2 ' 



(1) 



Notice that this metric is conformally equivalent to the flat 
metric on the unit disk i.e. ds^ 2 = £l 2 (x, 2/)^ s un j t di^- The 



boundary x 1 



y 



1 is a conformal boundary in that it 



is infinitely far away as measured in the hyperbolic metric. 
We will understand what role this space plays shortly. But 
before moving on, let us remember that we are interested in 
gapped topological liquids described by topological quantum 
field theories. We are perhaps not accustomed to thinking of 
these objects as conformal field theories, but while they have 
a much larger invariance group, they certainly belong to the 
more general class of conformal field theories. We will now 
build up to this result from a simpler result for the entangle- 
ment Hamiltonian of a half space. 

It has been known for some time, motivated by studies of 
Unruh radiation and Hawking radiation, that the spectrum of 
the reduced density matrix of a half space in a relativistic 
quantum field theory is related to that of the Lorentz boost 
generator that preserves the so-called Rindler wedge (see Ref. 
24). Focusing on the two dimensions x and t mixed by the 
boost, the Rindler wedge is the region of spacetime with 
x > 0, x < t, and x > —t. Considering the lightcone co- 
ordinates x + = t + x and x~ = t — x, boosts by velocity 

v in the x direction take x^ — > e ±x x ± (with e A = yjj^), 

and since the boundaries of the Rindler wedge are x + = 
and x~ = 0, we see immediately the boost only moves points 
around within the Rindler wedge. An alternate characteriza- 
tion of the Rindler wedge is simply the region of spacetime 
where the physics is totally controlled by the state of the half 
space x > at a fixed time, it is the "causal development" D 
of the half space. In fact, the reduced density matrix of the 



half space is nothing but a thermal state where the Hamilto- 
nian is the generator of Rindler boosts as shown in Ref. 24. 
This claim can be proved using a simple path integral argu- 
ment involving the Euclidean version of the Rindler boost as 
shown in Ref. 27, but it is already made plausible by noting 
that the map A — > A + 2ni leaves the coordinates unchanged. 
This is consistent with a quantum state that is thermal with 
respect to the generator of the boost. 

Returning to our original comments, this observation ac- 
counts for the radiation seen by accelerated observers, since 
their effective time evolution is generated by the Rindler 
boost. The Minkowski vacuum looks like a thermal state for 
the Rindler boost generator, and accelerated observers experi- 
ence such a state as an ordinary thermal bath with respect to 
their internal clock. The boost generator can be written ex- 
plicitly as K ~ xdt + td x , and we easily check that e XK x ± = 
e ±A i ± . The operator version of this generator acting on de- 
grees of freedom in the half space is K = J x>Q xH — VP 
where H and V are the energy and momentum density re- 
spectively. Evaluating this generator on the t = time slice 
provides an immediate connection between the entanglement 
Hamiltonian K and the physical energy density H. The only 
subtlety is the question of boundary conditions addressed be- 
low. 

Remarkably, the entanglement Hamiltonian is the genera- 
tor of a geometric flow in spacetime, and this flow may be 
interpreted as time evolution in Rindler space. The reduced 
density matrix of the half space is then a simple thermal state 
with respect to time evolution in Rindler space. If we change 
coordinates to x = a~ x e p cosh ai] and t = a~ 1 e p s'mharj, 
then the metric on the Rindler wedge takes the form ds 2 = 
e 2p (—dr] 2 + dp 2 ) and the curve p = has constant acceler- 
ation a. The spectrum of the entanglement Hamiltonian can 
then be found directly in the continuum by quantizing the low 
energy theory in this spacetime. A cutoff is necessary for 
large negative p (otherwise the effective temperature, given by 
1 / y/grjr), diverges as p — > — oo), a fact familiar from the study 
of black hole thermodynamics (the Rindler spacetime space- 
time approximates the near horizon limit of a black hole). 

The results for the Rindler wedge are already sufficient to 
prove essentially everything we want, but to work with com- 
pact subsystems we must introduce a little more technology. If 
we further restrict ourselves to a conformal field theory, then 
we have an additional result about the reduced density ma- 
trix of the d-ball as described in Refs. 25, 26. Those authors 
showed that the spectrum of the reduced density matrix of a 
disk is related to the spectrum of the conventional Hamiltonian 
of the conformal field theory but defined on the hyperbolic 
space H d x R. The proof is very similar to Rindler wedge re- 
sult, and actually follows from the result for the Rindler wedge 
by a conformal transformation. Ref. 26 contains details of the 
geometrical flow in the causal development D that is mapped 
to ordinary time evolution in the hyperbolic space. The pre- 
cise relation is that the spectrum of the disk density matrix is 
the same as the spectrum of a thermal state of the conformal 
field theory defined on the hyperbolic space. The radius of the 
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original disk sets the curvature scale of the hyperbolic space 
and the temperature of the thermal ensemble on that space. 

A very important check of this relation is that it correctly re- 
produces the divergences inherent in continuum entanglement 
entropy. In the hyperbolic setting, infinities arise because the 
hyperbolic space is non-compact. Hence the hyperbolic ther- 
mal entropy, which is the entanglement entropy, is infrared 
divergent due to the infinite volume. Cutting off the hyper- 
bolic space at large but finite size cures the divergence. As 
the explicit transformation shows, this large distance cutoff 
of the hyperbolic space is the same as the short distance cut- 
off in the field theory that renders the entanglement entropy 
finite. This is because the conformal boundary of the hyper- 
bolic space is mapped to the physical boundary of the orig- 
inal disk whose entanglement entropy we are interested in. 
The important message is that, although formally the hyper- 
bolic space is infinite, we must always provide a cutoff of that 
space at large distances. Thus we introduce a boundary into 
the space, and to compute the partition function properly, we 
must establish boundary conditions for the fields there. For 
example, boundary conditions must be considered to insure 
that the partition function is gauge invariant. 

A useful related construction allows us to obtain the entan- 
glement entropy of the disk (or 2-ball B 2 ) directly in terms of 
a partition function of the theory on a Euclidean sphere S* 3 . 
This result is also proven in Ref. 26. The proof proceeds as in 
the previous case by a judicious choice of coordinate transfor- 
mation and conformal transformation which maps the causal 
development of the disk into a spacetime whose Euclidean 
section is a sphere. Now the short distance divergence of the 
path integral on the sphere corresponds to the short distance 
singularities of the entanglement entropy. Since log Z gives 
the free energy, and since we will be interested in topological 
phases where the Hamiltonian is zero, the free energy directly 
determines the entropy (although in general one must be more 
careful as shown in Ref. 26). Before we proceed to give more 
intuition for these results, let us note that they have already 
been used to great effect in the holographic setting and in the 
study of subleading terms in the entanglement entropy of gap- 
less systems. These methods are also quite powerful in that 
they rigorously establish a boundary law for certain interact- 
ing gapless systems in more than one dimension. 

To develop an appreciation for these tools, let us give a few 
more details. We focus on the transformation to hyperbolic 
space since it is relevant for a finite subsystem, the disk. Also, 
we will describe the method in the context of 2+ 1 dimensions, 
but later we will use the 3 + 1 dimensional version as well 
(the result is completely general). Thus consider a disk, call 
it A, of radius R at fixed time t = in a conformal field 
theory in 2 + 1 dimensions. The reduced density matrix of the 
disk determines the physics within a spacetime domain called 
the causal development D of the disk. The boundary of this 
domain in spacetime has the form of two cones, one opening 
up and one opening down, that intersect at the disk at t = 
as shown in Fig. 1. These cones are light sheets that delimit 
the regions of spacetime reachable only by signals starting in 




FIG. 1: The disk A is colored blue. The rest of the system at fixed 
time is colored red and denoted B. The boundary causal develop- 
ment D is depicted with the wire outline cones; the solid cone is the 
forward boundary while the dotted cone is the past boundary. Every- 
thing inside the causal development D is determined by the state on 
A at time t = since no other region of spacetime can communicate 
with D. Figure adapted from Ref. 26. 

the disk at time t = 0. In terms of coordinates where the fiat 
space metric is ds 2 = —dt 2 + dr 2 + r 2 dfl 2 i _ 1 and the disk is 
r < R, the coordinate transformation which maps the causal 
development of the disk to the hyperbolic space is given by 

sinh (t/R) sinhu 
cosh u + cosh (t/R) cosh u + cosh (t/R) 

(2) 

In terms of these coordinates, the metric reads 

ds 2 = n 2 (-dr 2 + R 2 (du 2 + sinh 2 udfi^-i)) > ( 3 ) 

where the conformal factor fl 2 is (cosh it + cosh (t/R)) -1 . 
Now we simply remove the conformal factor via a confor- 
mal transformation , and we immediately have the confor- 
mal field theory in hyperbolic space. We can easily check 
that r — > R requires u — > oo. Using the same kinds of geo- 
metric arguments that we outlined for the Rindler wedge, one 
can show that the state of the conformal field theory on this 
spacetime is thermal under the Hamiltonian d T with inverse 
temperature 27ri?. 

To summarize, we introduced three powerful tools that pro- 
vide access to the reduced density matrix of special subsys- 
tems in terms of generators of geometric flows. First, the 
entanglement Hamiltonian of a half space in any relativistic 



Ref. 26 contains a minor typo: the vector C defining their conformal 
transformation should be (0, -1/(211), 0, 0). 
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FIG. 2: Visualizing the modular transformation that maps two copies 
of the solid torus B 2 x S 1 to the pair of interlocking solid tori on the 
right. Expanding the interlocking tori to fill all of space and gluing 
them along their mutual boundary produces S 3 . 



field theory is given by a certain boost generator. Second, 
the entanglement Hamiltonian of a ball in any conformal field 
theory is given by the generator of time translations in hyper- 
bolic space. Third, the entanglement entropy of a ball in any 
conformal field theory is given by the partition function of the 
Euclidean theory on S d+1 . Since we consider primarily topo- 
logical theories in this paper, all three results apply, but note 
that the half space result is especially general since it requires 
only Lorentz invariance. 



Warmup: topological entanglement entropy 

The simplest calculation that illustrates these tools is a com- 
putation of the entanglement entropy of a disk in topological 
liquids in 2 + 1 dimensions. As we already noted, topolog- 
ical field theories are certainly conformal field theories since 
they do not care about the metric of spacetime at all. To com- 
pute Z(S 3 ) we begin with the fact that Z(S 2 x S 1 ) = 1. 
This is the statement that a topological phase has a unique 
ground state on the sphere. The spacetime S 2 x S 1 may be 
cut open along the S 2 to yield two copies of the solid torus 
B 2 x S 1 . Since dB 2 x S 1 = S 1 x S\ Z(B 2 x S 1 ) = is 
a state in the Hilbert space of the torus generated by imag- 
inary time evolution. This state is normalized since 1 = 
Z{S 2 x S 1 ) = Instead of gluing the tori back to- 

gether directly, we first make an S modular transformation 
of one of the boundary tori which exchanges the two non- 
contractible surface loops. Now gluing the tori together yields 
S 3 as shown in Fig. 2. The modular transformation is imple- 
mented using the modular S'-matrix S£, and a simple calcu- 
lation gives Z(S 3 ) = (*|5|*) = 5§. Since S$ = 1/T> and 
using the fact that the ordinary Hamiltonian of a topological 
phase is zero, we conclude that S(B 2 ) = In Z(S 3 ) = — In V 
(the non-universal part has effectively been subtracted away) 
as originally shown in Refs. 5, 6. In fact, using the exponen- 
tially fast factorization of the density matrix for widely sepa- 
rated regions and the topological invariance of Stopo' we ma Y 
immediately establish that Stopo = — ng In I? where ng is the 
number of boundaries of the region of interest. 



BULK-EDGE CORRESPONDENCE 
Fractional quantum hall fluids 

Now we apply the technology introduced in the previous 
section the problem of the entanglement spectrum in frac- 
tional quantum Hall fluids. Consider the simplest class of 
topological fluids at filling fraction v = 1/m as described 
by Laughlin in Ref. 28. These states are described at low 
energy by an effective Chern-Simons theory for an emergent 
gauge field that can be used to compute ground state degen- 
eracy and quasiparticle statistics (see Ref. 29. We will now 
establish that the entanglement spectrum of a disk in this class 
of topological fluids is given by the thermal spectrum of a 1 + 1 
dimensional gapless edge mode. The argument is extremely 
simple. 

First, we map the reduced density matrix of the disk onto 
the thermal density matrix of the topological fluid living on 
hyperbolic space. To be completely clear, we recall that only 
the the space is hyperbolic, time is still an extra product di- 
mension. Second, if the radius of curvature is small, that is, 
if we asked about the entanglement entropy of a large disk, 
then the topological fluid does not care that it lives on curved 
space instead of flat space. However, there is one aspect of 
hyperbolic space that the fluid does care about: the existence 
of a (conformal) boundary. As far as the Chern-Simons theory 
is concerned, the fact that this boundary is infinitely far away 
according to the hyperbolic metric is irrelevant, and we may 
as well study the Chern-Simons theory on a disk. More prac- 
tically, we cut off the hyperbolic space at some fixed size and 
impose boundary conditions to regularize the path integral. 
As is well known, Chern-Simons theory on a manifold with 
boundary is not gauge invariant unless we add extra boundary 
degrees freedom [30]. Because the original electron model 
was gauge invariant, the reduced density matrix must also be, 
and thus we must include these edge degrees of freedom. But 
now we know the story very well, the edge degrees of freedom 
are those of a c = 1 gapless boson in 1 + 1 dimensions, and 
because the bulk is still fully gapped, these edge modes dom- 
inate the thermal physics. Hence the "low energy" part of the 
entanglement spectrum is simply given by the thermal spec- 
trum of the corresponding 1 + 1 dimensional conformal field 
theory. Of course, we can similarly conclude that the entan- 
glement spectrum of a half space in the Chern-Simons theory 
is given by the thermal spectrum of the same conformal field 
theory on an infinite line. 

Note also that we are free to cut off the hyperbolic space 
at very large size, much larger than the curvature scale, and 
this means that the spectrum of the edge theory can be arbi- 
trarily closely spaced since the spacing is set by the radius 
of the edge. For a particular microscopic model, we would 
have to cut off the hyperbolic space at a size set (via a confor- 
mal transformation) by the physical cutoff. Thus the details 
of the entanglement spectrum can depend on the microscopic 
cutoff and disk radius, but the existence of the gapless edge 
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FIG. 3: Left: physical edge spectrum of the model in Ref. 32 with 
a — 1 and R — oo. Right: spectrum of the approximate entangle- 
ment Hamiltonian with a = 1 and R = 20 (the spectrum is not too 
sensitive to the details of a and R and the overall scale is irrelevant). 
Note that we show both spectra on the same scale, but the spectrum 
of the entanglement Hamiltonian is spread over a much larger range 
to the linearly growing couplings. This also explains the difference 
in slopes for the edge modes as the edge at large r x has much larger 
couplings and hence a much larger speed. 

mode in the entanglement spectrum is robust. Finally, we ob- 
serve that the Renyi entropies are also easy to compute given 
the full spectrum. Conformal invariance fixes the form of the 
leading entanglement temperature dependence, and the sub- 
leading terms reproduce the findings of Ref. 31 that the topo- 
logical part of S a is independent of a. 

Topological band insulators 

These techniques can also be applied to topological insu- 
lators. Of course, for Chern insulators in 2 + 1 dimensions, 
the result is essentially identical to the quantum Hall case. As 
an illustration of the method, we directly calculate the spec- 
trum of a lattice version of the Rindler boost Hamiltonian. We 
take the lattice model in Ref. 32 whose low energy limit is 
massless Dirac fermions in 2d, and we approximate our ex- 
pression for the boost Hamiltonian by multiplying the hop- 
pings of that model with a linearly increasing function of the 
form a + r x /R where a is analogous to the cutoff we men- 
tioned above in Rindler space, and r x labels lattice sites in 
the x-direction. We cutoff this model at large r x to have a fi- 
nite Hamiltonian, so there is also an edge mode coming from 
the large r x boundary (this edge mode would be removed in 
a calculation of the spectra for the disk). The model remains 
translation invariant in the y-direction. We plot the spectra of 
this model in Fig. 3. 

Now consider the time reversal invariant Zi topological in- 
sulator in 3 + 1 dimensions with effective theory given by 
the 9 term 9 J F A F with 9 = 0, n. This is a rather simple 
field theory, but the 9 term is topological, so we may apply 
our procedure. We want to determine the entanglement spec- 
trum of a ball in this system, and as usual, we map it to the 
effective thermal problem in hyperbolic space. Cutting off the 



hyperbolic space at finite size once again leads to a physical 
boundary where we must specify boundary conditions in or- 
der to compute the partition function. There are two classes 
of boundary conditions we can consider. If we impose time 
reversal breaking (or particle number non-conserving) bound- 
ary conditions, then we will not in general have edge states. 
However, because we are asking about the reduced density 
matrix deep inside the bulk of a time reversal invariant system, 
the only sensible boundary conditions are those that preserve 
time reversal. 

We appeal to the robustness of the 9 = tt insulator to argue 
that if we do not break time reversal at the surface in our cut- 
off hyperbolic space, then we will have gapless edge states. 
Furthermore, these edge states will not be seriously perturbed 
by the fact that our system lives not in flat space but in hy- 
perbolic space as long as the curvature is smooth and the bulk 
gap persists. Examining the precise form of the entanglement 
Hamiltonian for the half space, we see that it describes a sys- 
tem where the effective gap increases from zero to infinity as 
a function of x, and so it is obvious that this system will have 
gapless edge states of the usual form in the non-interacting 
model. Microscopic interactions can be freely included so 
long as the system remains in the same phase and the low en- 
ergy effective theory is unchanged. Thus we have shown that 
the entanglement spectrum of a 3-ball B 3 in a Z-2 non-trivial 
insulator is necessarily gapless just like the thermal spectrum 
of a physical edge. 

A word of caution is appropriate here, since in the presence 
of interactions we cannot rule out the possibility of a surface 
phase transition. Of course, whatever the surface physics, the 
9 terms tells us that it must be gapless and must reproduce the 
1 /2 Hall response when time reversal is broken at the surface. 
Lorentz invariance permits us to seriously constrain the the- 
ory, but we can only really argue for surface Dirac cones in 
a weakly coupled description. If we have further information 
about the field theory, as in the standard non-interacting Dirac 
fermion model of topological insulators in 3 + 1 dimensions, 
then we can be quite precise about the half space entanglement 
Hamiltonian (we don't even have to formally send the gap to 
infinity to obtain a topological and hence conformal theory). 
On a final note, we could always consider strange boundary 
conditions, such as those that mimic the effects of disorder, 
but this would not be a good cutoff if the bulk preserved the 
relevant symmetries. These considerations do suggest an in- 
teresting question: if the bulk of a topological insulator is dis- 
ordered, could the entanglement spectrum track the change in 
the density of states at a physical surface due to this disorder? 

Fractional topological insulators in 3 + Id 

Precisely the same argument applies here for fractional 
topological insulators in d = 3 spatial dimensions. These in- 
sulators have a 9 term in their low energy effective action with 
fractional 9/ir as well as a topological BF theory (see Ref. 33 
for a brief definition) that describes ground state degeneracy. 
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In more familiar language, the low energy theory is decon- 
fined Z n gauge theory in 3 + 1 dimensions coupled to external 
electromagnetic fields via gapped fractionalized fermions that 
fill a topological band. Since these phases have protected edge 
states so long as time reversal is preserved, we may again ex- 
pect edge modes in the cutoff hyperbolic space or in Rindler 
space. These edge modes will dominate the low energy ther- 
mal spectrum and hence the universal part of the entanglement 
spectrum. The same caveats concerning the precise form of 
the edge modes applies here as well, but given a relativistic 
realization of the low energy effective theory, such as the one 
described in Refs. 34, 35, we can again be quite precise about 
the nature of the entanglement Hamiltonian. 

We also note that the sphere construction can be used to 
compute the topological entanglement entropy of a ball inside 
a fractionalized topological insulator. The topological entropy 
comes exclusively from the Z n gauge theory, and so we need 
to evaluate Zz n (S 4 ). If we normalize the partition function 
so that Z(S 3 x S 1 ) = 1 (one state on the 3-sphere), then 
we find that Z(S 4 ) = l/n. In fact, for any manifold M 
we have Z(M) = n bl( - Ar >~ 1 where b x (M) is the first Betti 
number of M, the number of non-trivial loops in M. This 
result follows, up to normalization, directly from the struc- 
ture of the Euclidean theory as the coupling is taken to zero. 
Using a lattice regulator, the flux through each plaquette is 
zero in the weak coupling limit, but there is a remaining free- 
dom to specify the flux through all the non-contractible loops. 
These configurations all occur with the same weight in the 
path integral and contribute a factor of n for each loop. Us- 
ing St p = In Z = — In n and setting n = 2, this formula 
reproduces the known result for 3 + 1 dimensional Z2 gauge 
theory obtained in Ref. 36. 

CONCLUSION 

We have established the bulk-edge correspondence for a 
wide variety of topological quantum fluids. This correspon- 
dence relates the spectrum of the reduced density matrix of 
a spatial subsystem in the bulk to the thermal spectrum of a 
physical edge. The entanglement cut becomes a physical cut. 
In addition to reproducing some old results in a unified frame- 
work, we have offered the first proof of the bulk-edge corre- 
spondence for fractional topological insulators. We do not be- 
lieve this result was seriously in doubt given the circumstantial 
evidence and intuition for the bulk-edge correspondence, but 
it is valuable to have a proof. Although we have addressed a 
wide variety of systems, we believe that our technique has not 
been exhausted. As an exact non-perturbative relationship be- 
tween the entanglement spectrum of simple subsystems and 
the thermal spectrum in simple spacetimes, we believe this 
technique has much to offer the study of entanglement in 
quantum many-body systems. As an example, it would be in- 
teresting to see if the ground state wavefunction can detect the 
full central charge of the edge or only the "topological cen- 
tral charge" given by c mod 8. Another subtle point is the 



issue of inversion symmetry described in Ref. 22 which may 
require a careful study of the boundary conditions or even the 
other Rindler wedge associated with x < 0. 

Although we only briefly mentioned it, our methods per- 
mit a calculation of the topological entanglement entropy of 
a ball in three dimensional topological phases. We have few 
physical candidates for such phases, and their entanglement 
properties are more complicated (because we must deal with 
flux lines, see Ref. 37 for a recent discussion), but we can reli- 
ably identify the subleading term in the entanglement entropy 
of a ball (and equivalent regions) from Z{S A ). 

Another interesting idea is that the mapping to hyperbolic 
space or Rindler space might be useful to numerically com- 
pute the entanglement spectrum of certain critical points. We 
would need to study a lattice model which realizes the con- 
tinuous quantum phase transition of interest on a lattice that 
mimics the appropriate geometry. We may be limited by our 
ability to numerically simulate such a model, but it should be 
possible for some models and would give direct access to the 
entanglement spectrum. On a related note, the mapping we 
have employed also suggests a way to compute entanglement 
entropies of critical theories very generally in terms of equiv- 
alent thermal problems. We suspect this can be used to make 
statements about the computational complexity of calculating 
the entanglement entropy for such systems on a quantum com- 
puter. Building on the work of Ref. 38, we are in the process 
of computing the entanglement entropy of balls at the 0(n) 
Wilson-Fisher fixed point. We end by noting that there are 
many other partitions of interest besides the spatial ones con- 
sidered here (e.g. Refs. 39, 40), so there is still a great deal 
to understand about the entanglement spectrum in quantum 
many-body systems. 
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